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1 Introduction 



A gauge/gravity duality model for SU(iV c ) gauge theory thermodynamics with a high tem- 
perature gluon plasma phase and a low T glueball gas phase and a first order transition in 
between has been presented in [H12]- The model contains a 5-dimensional metric ansatz and 
a scalar field with a bulk potential. A variant of the model has been developed in [3], in 
which the scalar field potential as a starting point is replaced by the beta function of the 
boundary field theory. 

The purpose of this article is to apply the method in [3] to a boundary field theory with 
the model beta function [U [5] 

/3(A) = -/3 A 2 f 1 -^)' A = N c9 2 , (1) 

where /?o > is a parameter and A = A* is an infrared (IR) fixed point. What is the finite 
temperature phase structure of this theory? At high temperatures one expects the existence 
of a "quark-gluon plasma" phasa^j conformal at T — > oo. Due to the IR fixed point one 
knows that there are no massive glueball states in the theory. Thus there is no low T glueball 
phase but, instead, a conformal "unparticle" phase of massless particles. We shall show that, 
indeed, this picture follows naturally from the model in [HE]. Furthermore, the transition is 
of first order for /3oA* > 6.58 and a crossover for smaller values. 

More formally, at T = and T = oo we have two pure AdSs spaces with radii £ir and 
£uv and the transition we study is driven by varying T. As we shall see, the parameter /3o 
is essentially the ratio of the radii: 

A,A, = §log^. (2) 

Beta functions of the type in Eq. ([!]) but somewhat more complicated have recently become 
a central topic of research in connection of technicolor theories [SI El El E] and lattice Monte 
Carlo studies thereof pH HH Q21 EG9 O EES QI] . The results in this article should set the 
stage for the study of the thermodynamics of these more general theories. Is it possible to 
use lattice Monte Carlo to search for the two-phase structure predicted here and to determine 
the order of the transition? 



2 The gravity sector of the model 



We summarize first the gravity sector equations of the model in [TJ [3] . The model starts from 
a metric ansatz 

ds 2 = b 2 (z) -f( z )dt 2 + d* 2 + ^- (3) 

J\ z ). 

plus a scalar field (fi(z) = log A(z), where the four functions b(z), f{z) in the metric, the scalar 
field 4>(z) and the potential V(4>(z)) in the gravity action (in standard notation) 



(4) 



1 Realistic theories of the type have fermions and we assume that we can apply the model in [3] to 
them. 
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are determined from the four equations (b = b'(z), etc.) 

6 P + 3 6 + 3 67 = 7 

b 2 b , - 9 

4 + 3- = 0, 
/ 6 

/3(A) = 6^, X(z) = e*M ~ 5 2 iV c 
First, from ([8]) it follows that, for the beta function ([I]), 
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1 i A 

- - log - 
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where we have scaled (hats will be removed after Eq. (|16|)) 



Q + ilog(/3 Q-l), 
Po 



A — t— , /3o — A*/3o 

A* 



and where we have introduced [3] 



Q = 



1 



Further, to solve the second equation one introduces 

w = -b/b 2 , 

and solves 

W(X) = W(0)exp(-lJ Q X dX^)=W(0)e 



|/3o(A-|A 2 ) 
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(11) 
(12) 



(13) 



Given W one can immediately write down [U [5] the scalar field potential V((f>) corresponding 
to f(z) = 1 in ([5]) (from now on hats are removed): 



V(X = e+) = ^- 2 e> e-P°^ 2 [l - i f3 2 X 2 (l - A) 2 ] . 



(14) 



In the UV (A = 0) and IR (A = 1) limits we have two AdSs spaces with radii £ = £tjv 
and £fr = £uve _2/3 °^ 9 - Eq. ([2]), giving the parameter /3q as essentially their ratio, follows 
immediately from this. 

When / is the solution of Q, (jSJ) gives the potential 



V(A) = !2fW 2 



3A 



f 
b 



(15) 
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From Eq. (12) one further obtains, defining the energy scale A = bo/C, 

dQ 



Az 



q (A)Q - iy'^ +i 



foQe Q 9 <? + W^. 



Finally, from Eq. we solve 



m = i 



dz 

o W 



I 



dz 



(16) 



(17) 



where two constants of integration are fixed by /(0) = 1 and /(%) = 0. 
Thermodynamics now is obtained from the key relations 
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For the temperature we thus have 



dz 



4vrT 



b J > 
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(18) 



(19) 



or 



A 
4^rT 



= 3Q 



(A<3-1) 4 '« 

The derivative dT/dQ will be very important in what follows and this is computed to be 

dT (3 Q /4vrT 



(20) 



TdQ (3 Q - 1 \ bW 
Full thermodynamics is then given by 
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where 6(Q) is given by @ and T(Q) by (|20j). 
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3 Analytic approximations in the ultraviolet and infrared 

Final results will have to be computed numerically, but it is very useful to have analytic 
approximations in the ends of the phase space, UV (Q — > oo) and IR (Q — > 1//3q)- 
The leading UV terms are very simple: 



Oo Az oqz A 



(27) 



further corrections can be easily evaluated. The UV limit of bulk thermodynamic quantities 
is easiest to obtain from 



£ 3 (o(Q^oo)/6 ) 



■ 7T 



T 3 4G 5 (T(Q oo)/A) 3 4G 5 
since the UV limit here is nothing but 1/(z/jT) 3 = 7r 3 . Equivalently, 



V_ 

rp4 



£ 3 vr 3 



4G, 4 ' 
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Normalising to ideal gas pressure p = <? e ff7r 2 r 4 /90 one has 

£ 3 2 



4G, 



457T 



5cff- 



(28) 



(29) 



(30) 



The conformal IR region is the region near /3$Q — 1 = where the powerlike behavior of 
b(Q) drives it to zero. For b and W one can trivially expand: 



b 
W 



b e?o ((3 Q-l)Po 



Po 



D + 
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(/3 Q-1) 2 + ... 
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(/3bQ - l) 2 + O((3 Q - 1) 



(31) 
(32) 



For Az and T/A the situation is more subtle: to ensure convergence of the integrals (16) and 
( 20 ) at the lower limit one must partially integrate in the power of /3qQ — 1 to decrease the 

(33) 



exponent. One partial integration gives the leading power. For Az one has 

Az 

where 



e ^°(A)Q-i) *+/*(#,)][! + 0(A>Q-l)], 



dQ 



d 



(34) 



l/fh (A)Q - l) 1 ^ dQ 

is an integral arising from partial integration. This converges for /3q > 1. If ft < 1 more 
partial integrations are needed to obtain finite subleading terms. Anyway the leading term 
(33) implies that the fields always extend to z = oo in the IR limit. 



For T the leading term is 
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-e^ 9 (ftQ 
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(35) 
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where 



o) 



dQ 



/3 Qe 



-4Q- 



9Q 



A) >4, 



(36) 



converges if (3q > 4. Using this in (26) one obtains for the sound velocity in the conformal 
region 



1 



l_ 4e/3o + 9 I T (p )(p Q Q - l)W 



(37) 



For the pressure the next-to-leading order expansion in the IR conformal domain is (see 
Figs. [3] and [§ 



P C 3 7T 3 jgfl 

— j = e 3 

T 4 4G 5 4 



4 , 20 



l + | e /3o- 9 J T (/3 )(^oQ-l)^ +O(/3 Q-l) 



(38) 



The overall magnitude here is the same as at large T but suppressed by the factor e" 



4 The two phases 

To clarify the structure of the solution, the behavior of T{Q) for a few values of (3q is shown 
in Fig. |l[ For any value of /3q, T vanishes T(Q)/A ~ (/3qQ — l) 1 ^ at the infrared endpoint, 



see Eq. (g. For Q > l/p the limit is T(Q)/A ~ e Q , Eq. (27|. In between, /3 = 6.58 
separates two qualitatively different behaviors. For (3q < 6.58 T(Q) grows monotonically 
with Q, dT/dQ > always. For (3q > 6.58, T(Q) has a maximum and minimum, i.e., two 
zeroes of dT/dQ and a range where dT/dQ < 0, (? s < so that the system is mechanically 
unstable. 

To see concretely what this implies for physics, consider separately the two cases /3o = 8 
with two zeroes of dT/dQ and /3q = 2 with dT/dQ > 0. 
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Figure 2: (Left panel) The pressure for /3 = 8. Note the positive peak at Q = 0.137, caused 
via Eq. (23) by the increase of T(Q) starting from the IR limit 1//3q = 0.125. (Right panel) Field 



configurations in the conformal region for /3o = 8. Note that the sound velocity vanishes at the same 
point Q — Q max — 0.137 where p'(Q) and T'(Q) vanish. 
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Figure 3: Plot of p/T 4 vs T/A in the critical region for a non-monotonic T(Q). There is a stable high 
T gluon plasma phase with a metastable supercooled (SC) branch and a stable low T unparticle phase 
with a metastable superheated (SH) branch. The metastable branches terminate at the points where 
sound velocity squared becomes negative; between these points the system is mechanically unstable. 
There is a first order phase transition at T c = 1.084A. The unparticle phase extends to T — to the 
value 7r 3 e _16 / 3 /4 = 0.037, the gluon plasma phase at T ^> T c to 7r 3 /4 = 7.75. For normalisation to 
ideal gas, multiply by the factor in Eq. (30). 



Fig. [2] shows for (5q = 8 a plot of p(Q) and a blow-up of the far infrared region. The two 
zeroes of T'{Q) correspond to the two zeroes of p'(Q), a maximum and a minimum of p(Q). 
The physics is seen from the plot of pressure vs. T in Fig. [3j The correspondence between 
Figs. [2] and [3] is as follows: 

• The range 1/Po < Q up to the first maximum Qmax of p(Q), at which c 2 s = 0, corre- 
sponds to the unparticle phase in Fig. [3j Only the part of this up to Q c \ is thermody- 
namically stable, the rest is a metastable superheated branch. 
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The interval Q meiX < Q < Q mm in which p(Q) decreases in Fig. [2] corresponds to the 
unstable phase with c 2 s < in Fig. [3} 

The interval of Q > Q mm in which p(Q) again increases in Fig. [2] corresponds to the 
quark-gluon plasma phase; only the part of this above Q C 2 is thermodynamically stable, 
the rest is a metastable supercooled branch. 

The phase transition takes place at T c = 1.084 A. The latent heat is numerically 

Ae £ 3 

4.55, (39) 



T c 4 4G 5 



where (30 ) gives the normalisation. A fit to nearby values is L/T 4 = 3.86(/3q — 6.58) ' 4 ' 2 ' 



If one now decreases /3o, the first order transition in Fig. [3] becomes weaker and terminates 
at Pq = 6.582 at T = 0.991 A and p/T 4 = 0.134 (which decreases to 0.0963 at T = for this 
/?())• Below that the transition is continuous. The situation for @q = 2 is shown in Fig. |4| Since 
dT/dQ > there is no range with negative c 2 s and the unstable phase in Fig. [3] disappears. 
Related to this, there is no discontinuity in p'(T), the two phases are continuously connected. 
Fig. [4] also shows the interaction measure 

e — 3p m d p , , , 

(40) 



T 4 dT T 4 ' 

Its maximum at T = 1.81 A could be used as an estimate of the crossover temperature at this 

A)- 




Figure 4: Plot of p/T 4 and of the interaction measure (e— 3p)/T 4 (multiplied by 5) vs T/A for /3 = 2 
with a monotonic T(Q). The two phases are continuously connected. The unparticle phase extends 
to T = to the value tt 3 e" 4 / 3 /4 = 2.043, the gluon plasma phase at T > A to vr 3 /4 = 7.75. For 
normalisation to ideal gas, multiply by the factor in Eq. (30). 
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5 Conclusions 



We have, in this article, studied the thermodynamics of a field theory with the beta function 
([!]) with an infrared fixed point using a gauge/gravity duality model [HE1E]. In the two limits 
T = and T = oo we have conformal theories with pure AdSs spaces as gravity duals. The 
crucial parameter /?o of this model is essentially the ratio of the radii of these spaces, Eq. Q. 
Varying T drives a phase transition which in this model is of first order for £uv/£ir > 4.32, 
continuous otherwise. 

The computations here used Eqs.Q-Q to determine the metric functions, scalar field, 
and the phase structure [3]. However, we have also verified that the same phase structure is 
obtained by integrating numerically the Einstein equations ([5])- ([7]) with the scalar potential in 
(14) and by computing the equation of state as in [2]. Quantitative details change somewhat, 
for example, the critical value of /3oA* is 4.85 instead of 6.58 above. 

It has been a very challenging task to identify four dimensional gauge field theories with 
beta functions of the type Q. These need fermions in special representations [6] and the 
proof requires extensive numerical effort [TUJ [TTJ [J31 HH [JSJ [TB] . It is clear that the finite 
temperature properties of these theories will be studied but the path to clarification of the 
phase structure will be long. 

Theoretically it will be interesting to study the finite T phase structure of theories with 
a walking technicolor beta function, for which the IR fixed point is only approached. The 
theory then is confining and there is no conformal unparticle phase at T = 0. The model in 
this article can be used for that purpose. 
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